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Exciton-Polariton Quantum Gates Based on Continuous Variables 
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We propose a continuous variable analog of quantum controlled-NOT gates based on a system 
of exciton-polaritons in semiconductor microcavities. This can be realized by the engineering of 
parametric interaction between control and target polariton modes, which can be varied in time. 
As an explicit setup we use a system of dipolaritons, which allows for enhancement of parametric 
interaction by auxiliary classical fields and scalable multigate system realization. The calculated 
fidelity is shown to exceed 99% for realistic system parameters. 
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Introduction .—Quantum controlled NOT (CNOT) 
logic gates are universal elements in quantum com¬ 
putation, in principle allowing the implementation of 
any quantum algorithm (when supplemented with sin¬ 
gle qubit rotations). Their construction in physical sys¬ 
tems is challenging since one requires a system with both 
limited dephasing and strong nonlinearity. Early realiza¬ 
tions of CNOT gates made use of nuclear magnetic reso¬ 
nance in molecules [ij or post-selection in linear optical 
systems 0 , which were later reduced in size on photonic 
chips [^. More recently, demonstrations of quantum 
gates and small quantum circuits were achieved using 
a variety of systems, including ion traps , nitrogen 
vacancy centers 0, and superconducting qubits (ol-fl^. 

Semiconductor systems have long been valued in in¬ 
formation processing for their compact sizes, which 
are particularly important when one aims at building 
circuits with large numbers of gates. In particular, 
semiconductor microcavities containing quantum wells 
are quasi-two-dimensional structures with micron thick¬ 
nesses. These structures also offer the opportunity to 
hybridize the properties of photons and excitons, gen¬ 
erating new exciton-polariton quasiparticles with decay 
time exceeding tens of picoseconds. The presence of an 
excitonic component facilitates nonlinear interaction be¬ 
tween polaritons, which has led to experimental observa¬ 
tion of quantum optical effects such as squeezing [13| and 
non-classical correlations 1^ for a coherently driven po¬ 
lariton system. Additionally, macroscopically populated 
polaritonic modes were suggested to mimic the two-level 
qubit system and analog CNOT gate 0. At the 
same time, currently existing polaritonic samples do not 
possess the strong single polariton nonlinearity required 
for the conventional blockade mechanism 17|, and most 
quantum effects at the weak single polariton nonlinearity 
level need to exploit interference effects for realization of 
an unconventional polariton blockade 00 or genera- 
tion of entangled states [2^ . 

The aim of this paper is to introduce an alternative 
route towards construction of polariton quantum CNOT 
gates. To avoid the need of single particle control and 
detection, we choose to encode information in the con¬ 


tinuous amplitude and phase variables of polariton fields. 
First, a continuous variable quantum CNOT gate opera¬ 
tor is realized by engineering a two-mode parametric in¬ 
teraction Hamiltonian, which can be controlled in time. 
Second, it is well-known that it is difficult to observe 
nonlinear effects when small numbers of polaritons are 
involved. For this reason we exploit a mechanism of non¬ 
linear enhancement that uses a macroscopic density to 
amplify scattering processes, while quantum information 
is maintained in separated low density quantum modes. 
The proposal allows to achieve CNOT operation fidelity 
as high as 99% for realistic system parameters. 

Finally, we consider the scalability of our scheme, 
where multiple quantum gates can be cascaded one af¬ 
ter another, in principle allowing the construction of ar¬ 
bitrary algorithms. Typically one imagines a quantum 
circuit as a network of quantum logic gates separated in 
space [i^. This requires that signals travel spatially be¬ 
tween distant nodes of the network. However, spatially 
propagating polaritons would experience losses as they 
scatter with disorder and experience dispersion. While 
superfluidity [ii-0 and soliton wavepackets [1^ have 
been shown to overcome these effects in microcavities, 
they imply a macroscopic classical polariton state, which 
cannot itself encode quantum information. In the present 
proposal we use a reciprocal space encoding of quantum 
nodes, negating the need for any spatial propagation. 
Polariton modes are distinguished by different momenta 
and logic gates exploit momentum conservation rules to 
connect particular modes. Rather than being physically 
fixed, the gates are enacted by the application of a known 
pulse sequence, which could be controlled by a spatial 
light modulator. This brings the additional feature of 
being able to reconfigure the quantum circuit. 

Definitions. —The theory of quantum information with 
continuous variables [Hj is well-developed and uses the 
same fundamental features of quantum computation, 
namely superposition and entanglement, to achieve sim¬ 
ilar aims to qubit based methods. Working with contin¬ 
uous variables, the analogous definition of the quantum 
CNOT gate is given by the operator CN = 0, 

which acts simultaneously on two quantum helds, di and 
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0 , 2 , where the amplitude and phase operators are defined 
by qn = {an + al) /v^ and = -i (a„ - aj,) /V2. The 
CNOT gate can also be defined by its action on ampli¬ 
tude eigenstates CN|gi,q 2 ) = \qi,q 2 + qi), demonstrat¬ 
ing analogy with CNOT gates: the first quantum field 
Si acts as a control field, which adds to the amplitude of 
the target field 0,2 ■ 

Single gate scheme .— To realize a single CNOT gate, 
let us consider the Hamiltonian: 

a = a ('0*^0102 -I- 0^S|S2) — J (a|S2 -|- SJSi) . (1) 


It represents a parametric scattering from macroscopi- 
cally occupied modes 0 into two quantum modes Si and 
S 2 , which have an additional linear coupling between 
them. If we imagine that Si and S 2 represent polariton 
modes with equal and opposite wavevectors, then the first 
term in the above Hamiltonian can be arranged making 
use of interbranch scattering schemes such as those stud¬ 
ied theoretically in standard planar microcavities 28| and 
parametrically driven cavities [ 2 ^, or experimentally in 
triple microcavities [s^- In such cases the field 0 would 
be a polariton field with zero in-plane wavevector, and 
a is a constant parameterizing the strength of polariton- 
polariton interactions. We will assume that 0 can be 
controlled optically, and ideally the linear coup ling term 
J could also be controlled. Both resonant [SlJ and non¬ 
resonant excitation have been used to optically gen¬ 
erate polariton potentials, lea ding to the optical control 
of polaritons in real-space and creating the de¬ 

sired linear coupling J between the modes di and 0,2 ■ 
Since the amplitude and phase of 0 can be tuned, let 
us consider the case a0^ = — J. The phase reference of 
the second mode can be changed 0,2 m 2 following the 
phase of an input. Then, the Hamiltonian becomes. 


a = J (^—iaia2 + iajdj — ia\a2 + mjdi^ = 2JqiP2, 

~ ( 2 ) 

and the corresponding evolution operator reads U = 
By switching on and off the cou¬ 
plings the evolution time could be set to t —>■ tq = h/2J, 
in which case the unitary evolution operator corresponds 
to the quantum CNOT gate. The result of CNOT gate 
operation can then be monitored by measuring the am¬ 
plitude of a target field [s^. We note that while the 
interaction energy between a pair of polaritons 2 q; may 
be limited, it is the quantity q; 0^ that determines the 
relevant coupling strength [37|, where weak interaction 
is effectively amplified by the macroscopic field 0. 

Gate fidelity .—In the ideal case of unitary evolution 
with precisely timed parametric interaction, which we 
described before, the polaritonic system can work as a 
perfect continuous-variable CNOT gate. However, in real 
systems we identify the main mechanisms responsible for 
imperfect gate operation to be the decay of polaritons 
and the imprecision of the interaction constant control. 


In the presence of decoherence, the evolution of any 
expectation value is given by: 


irm. 

dt 


= A,n 






(3) 


where T denotes a dissipation rate of polaritonic modes, 
and the Lindblad superoperator is defined as = 

2a^Aa — d^dA — Ha^a, and Lindbladian £^ta[^] = 
2d^dHd'l'd — d^dd^dA — Ha^aa^a corresponds to a pure 
dephasing with rate Tp. For the sake of generality, we 
also introduced an additional incoherent pumping at rate 
P with conjugate Lindbladian which is respon¬ 

sible for an incoherent replenishing of polaritonic states, 
and may arise from a presence of a thermal reservoir. 

Using Eq. we derive a closed set of evolution equa¬ 
tions for the amplitude and phase expectation values: 


-( 92 ), 


i(p-r-rp),., 

=-2- 

■.dih) o-r/^ N , *(^-r-rp, 

^- 

^d(pi) _ T/~x i (P — ^ ^ p) ,, 

= -2iJ(p 2) H-^-^(Pi), 

.^d{p2) *(P-r-rp),., 

=-2- 


(4) 

(5) 

( 6 ) 
(7) 


which are readily solved analytically. One sees that for 
P = r the mean-field amplitudes of the polaritons no 
longer decay, while in the regime P > F one should ac¬ 
count also for nonlinear losses 3^ to prevent the ampli¬ 


tudes growing indefinitely. Furthermore, one can write 
a closed set of evolution equations for the second order 
correlators {{q\), ($ 192 ), ( 92)1 (Pi)i etc.). These equa¬ 
tions are given in the Supplemental Material [s^, to¬ 
gether with their analytical solution. 

To assess the performance of the quantum gate, 
we consider a set of displaced squeezed vacuum 
states as inputs Pm = |0i)(0i| ® 102 ) ( 02 |, with 

I0n) = i)(C,p“)5(r) |0), where = 

exp {{a^n)d\ — (a“)*a„) is the displacement operator 
providing the amplitude q™ and phase of the in¬ 
put mean-fields (a“). S{r) = exp is the 

squeezing operator with squeezing parameter r, and 
|0) is the vacuum state. The ideal output state is 
Pideai = 1 9 “,qT + qT) ( 9 “^+ 'J'™!- The actual state 
p{t) obtained in the presence of dissipation and incoher¬ 
ent pumping is characterized by its fidelity with Pideai: 


F {pidea,h p{t)) = Tr 




( 8 ) 


The fidelity can be calculated from the covariance ma¬ 
trix [d^, li^, which is directly obtained by knowing all 






















3 



FIG. 1: (color online) Fidelity of a quantum CNOT gate for 
increasing dissipation rate. The different colors correspond to 
different amounts of input squeezing (magenta: r=0, green: 
r=0.5, and blue: r=l). The solid curves show the range of 
fidelities for total intensity increasing from (a:}^ = 0 to {x)^ = 
10, where the highest fidelity is obtained for the squeezed 
vacuum input. Dashed curves show the same result in the 
presence of an additional incoherent pumping, chosen with 
rate P = F. b) Dependence of the minimum fidelity of the 
input squeezing r and total intensity {x)^. The different sur¬ 
faces correspond to different values of F/J = (0.01, 0.02, 0.05). 
Smaller values of F/ J leading to very high fidelities are ex¬ 
pected in samples with high quality factor. 


the second order correlators. The fidelity of the polari- 
tonic CNOT gate as a function of decay rate is shown 
in Fig. [IKa), for different values of the squeezing pa¬ 
rameter r. For each value of r we minimize the fidelity 
over input states with a fixed maximum total intensity 
= (qi)^ -I- (pi)^ -I- (92)^ + (P2)^. We assume the op¬ 
timum operation time tq = h/2J. The fidelity is seen 
to drop from unity (corresponding to a perfect CNOT 
gate) monotonically as the ratio of the decay rate to the 
coupling strength F/ J is increased, or when input states 
with a higher squeezing or total intensity are considered, 
as illustrated in Fig. [Hb). 

The decrease of the fidelity when operating with more 
highly squeezed states is very natural. In phase space, 
the Wigner function of the unsqueezed state is a circular 
Gaussian, while the infinitely squeezed state is a thin line. 
The decay has the effect of smearing any squeezed state 
into the unsqueezed state with (g^) = (p^) = 1/2. In fact 
the state is driven toward the unsqueezed state at a rate 
proportional to 1 /2— (g^) (as can be seen from writing the 
correlator evolution explicitly), that is, highly squeezed 
states are most quickly deformed in phase space. 

The dashed curves show the fidelities calculated with 
a non-zero incoherent pump, which attempts to compen¬ 
sate the losses in the system (P = F). Interestingly, in 
this case the fidelity no longer depends on the intensity 
of the initial state, however it becomes worse due to the 
incoherent pump. While the incoherent pump can com¬ 
pensate fully the loss of mean-field amplitudes, it can not 
compensate the loss of quantum correlations caused by 
dissipation. Allowing for different values of P, one finds 
that the optimum fidelity appears for P = 0. 



Fp/J 


Tf/Xo 


FIG. 2: (color online) Range of GNOT fidelities calculated 
for input states of total intensity increasing from (a:)^ = 1 to 
{x)^ = 10, accounting for a pure dephasing and non-ideal gate 
pulse, a) Fidelity as a function of pure dephsing rate for input 
states with squeezing parameter r = 0, 0.5, 1. b) Gate fidelity 
in square J{t) pulse operation mode plotted as a function of 
signal edge time to optimal operation time ratio, xf jto. 


Next, we plot the dependence of gate fidelity on the 
pure dephasing rate. The results are shown in Fig. [2Ka) 
assuming the decay rate F/J = 0.02, input states with 
r = 0,0.5,1, and varying total intensity ranging from 
(x)^ = 1 to (a;)^ = 10. These suggest that a low level of 
pure dephasing Fp/J < 0.05 is required for high quality 
CNOT operation. 

Additionally, we study the influence of imperfect tim¬ 
ing of the gate pulse on the fidelity of the continuous 
variable CNOT gate. We start considering a square pulse 
with finite pulse edges described by J{t) = Jof{t)[l — 
f{t — To)], where f{t) = [exp(—t/ry) -|- 1]“^, with Tf be¬ 
ing pulse edge time, and Jq denotes time independent 
interaction constant. The fidelity as a function of pulse 
edge time is shown in Fig. [21(b) for an ideal case of negli¬ 
gible losses. The results suggest that an accurate control 
of the gate pulse is required at a time scale below h/Jo, 
where pulses with sharp edges, t//to < 0.3 do not con¬ 
tribute to degradation of the gate fidelity. 

Finally, we discuss the experimental feasibility of the 
proposed scheme, with the main parameters being the 
polariton lifetime, nonlinear coupling constant, pure de¬ 
phasing, and characteristic switching time. Polariton life¬ 
times Tdec in the range of hundreds of picoseconds have 
been reported recently [3 j with the corresponding 
decay rate being in the fieV range. When controlled with 
a classical field, nonlinearities can be tuned to be in the 
sub-meV range. Taking J = 0.2 meV and Tdec = 65 ps 
gives the ratio of F/J « 0.02, and sets the characteris¬ 
tic switching time to a picosecond range, which shall be 
achievable with an optical control of the coupling con¬ 
stant. The pure dephasing in polaritonic system was 
estimated to be at Fp = 0.2 /reV level [^, yielding 
Fp/J = 10“^. This makes high fidelities of 99% feasible 
for relevant input states of {x^) < 10 and r « 0.5. 

Multi-gate system: dipolariton setup. —We now discuss 
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FIG. 3: (color online) Scheme for a CNOT gate implemented 
between a selected pair of modes. The plot shows 2D disper¬ 
sion of three polaritonic modes, each split in two branches for 
different polarization components. 


a scalable scheme of CNOT gates, which requires a set of 
multiple quantum modes and the ability to apply succes¬ 
sive gates between chosen pairs of modes. Toward this 
aim let us consider a dipolariton system , in which two 
types of exciton modes (direct and indirect) are coupled 
to a cavity mode in a microcavity resulting in three dis¬ 
persion branches as illustrated in Fig. [3l The advantage 
of the dipolariton system is the freedom in varying the 
mode energies via an applied electric field, although sim¬ 
ilar setups could be imagined in triple microcavities . 

Each dispersion mode exhibits a linear polariza¬ 
tion splitting between the transverse-electric (TE) and 
transverse-magnetic (TM) polarizations [i^. Let us take 
the TM polarized modes of the middle branch as the rel¬ 
evant quantum modes of our system. A degenerate set 
of these modes can be distinguished by different orienta¬ 
tions of in-plane wavevector and are illustrated in green 
in Eig. [3l In the following we will choose to work in the 
frame rotating at the energy of this set of modes. 


The Hamiltonian of the multimode system up to an 


energy shift reads (see details in [36|, sec. 3): 


n = AV'up.-i-l/'LP.+a2 +^i + + -SV'lp -'0LP,-a2 -«i,- 

’ (9) 

-I- C'V'up,+V'LP,-a 2 ,-i-®i.- + ^V'lp.-V'lp-S l+Si,-!- + h.c. 


where A, B, C, and D are effective interaction constants 
which depend on Hopfield coefficients and bare interac¬ 
tions. 

The quantum modes can be re-written in terms of their 
TE and TM polarized components using the transforma¬ 


( ^n, + ^n,— ) — 


1 / 


o2i0n 


^2 V 


( Sri. TMUra.TE ) , 

( 10 ) 

where (f>n denotes the angle of the mode n in reciprocal 
space. 

To arrange for a CNOT gate between an arbitrary pair 
of modes, oi^tm and 02 ,tm, we consider a quite specific 
but fully feasible macroscopic excitation of classical fields 
in the lower and upper polariton branches, as illustrated 
in red in Fig. [3l In the upper branch, we choose a cir¬ 
cularly polarized field '0up,-i- at the same wavevector as 
Oi.TM- In the lower branch, we choose a cross-circularly 
polarized field V’lp,- with the same wavevector as 
together with a linear polarization (characterized by a 
superposition of both circular components ip'LP,+ and 
iphP.-) at the same wavevector as 02 ,tm- Since we have 
chosen classical fields with the same in-plane wavevec- 
tors as oi^tm or 02 ,tm, momentum conservation rules al¬ 
low only a coupling of the chosen quantum modes by the 
classical helds. In particular, two types of scattering pro¬ 
cesses appear, where interactions between parallel spins 
are considered. First, the (T+ polarized component of the 
pumping provides an effective two-mode squeezing inter¬ 
action, allowing for the scattering processes indicated by 
the solid arrows in Fig. [3] (the reverse scattering process 
also occurs). Second, the cr_ polarized component of the 
pumping results in a process where the quantum modes 
are exchanged, by transferring their momentum to the 
classical field. This process is indicated by the dashed ar¬ 
rows in Fig. [3 (and again the reverse process also occurs). 
Finally, the classical field energies can be chosen such 
that the relevant scattering processes are resonant only 
with the TM polarized states. Neglecting off-resonant in¬ 
teractions with TE states (see discussion in [36|, sec. 4) 
gives the Hamiltonian: 


H = iJa\ r,,,, -I- iJal 




where J can be made real by correct choice of the am¬ 
plitudes and phases of the classical fields ([1^, sec. 3), 
allowing for CNOT operations acting on eight different 
continuous variable modes. 

Conclusion .—We presented a scheme for quantum 
logic gates based on exciton-polaritons in semiconductor 
microcavities. Unlike previous schemes, we operate with 
continuous variables that avoid the necessity of operat¬ 
ing with a definite number of polaritons, and make use 
of an effective amplihcation of nonlinearity in the system 
based on the coupling of quantum modes with macro- 
scopically occupied classical states. Using these ingre¬ 
dients, a quantum optical treatment of decay processes 
predicts fidelities in excess of 99% for existing microcav¬ 
ities. Additionally, we proposed a way for a construction 
of scalable networks of polaritonic gates. The experimen¬ 
tal demonstration of our proposal would not be reliant on 
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single-photon detection and would require standard ho¬ 
modyne detection measurements (this would depend also 
on the implementation of future error correction proto¬ 
cols). 
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1. Possible implementation of an experimental 
detection scheme 

The characterization of a quantum optical process can 
be achieved by measuring its effect on input coherent 
states [I|. Gaussian coherent states are fully character¬ 
ized by their covariance matrix, which can be experimen¬ 
tally accessed using homodyne detection Q. An impor¬ 
tant ingredient of such a technique is the availability of a 
classical local oscillator, with identical frequency to the 
measured modes. While this is not immediately available 
in our system, a modified detection scheme can be im¬ 
plemented using an additional interference between the 
lasers driving upper and lower branch polaritons. 

For the described system, the relevant two-mode co- 
variance matrix is given by 

Vij = i {xiXj + XjX,) - (xi) (xj ), (SI) 
where q 2 ,P 2 )- 

From the interference of the modes di and 0,2 on a 
beam-splitter and the application of phase delays, one 
has access to the fields 03 = (Si -I- 0 , 2 ) 0,4 = (Si — 

S 2 )/'\/ 2 , S 5 = (fSi -I- S 2 )/v^, Se = (fSi - 0 . 2 )!'/^ [see 
Fig.IHUa)]. 

The fields Si all oscillate at the same frequency, which 
is midway between the frequency of the laser driving the 
upper polariton branch and the frequency of the laser 
driving the lower polariton branch. Let us consider in¬ 
terfering any of the fields Si on a beam-splitter with a 
local oscillator field of the form qlo -I- e“**^*). This 
form of local oscillator can be attained from the superpo¬ 
sition of the lasers driving the lower and upper branches, 
where ano is the (complex) field amplitude and is the 
frequency difference with the middle polariton inodes. 
The photocurrents obtained at each of the output ports 
are [see Fig. ISlT bi]: 

= 2 + '^LO 

X (S, + aLo(e*''‘ + e-*''‘)), (S2) 

n_, = i(st-<o (e*^‘+e-“)) 

X (S,-aLo(e*''‘ + e-*''‘)). (S3) 

The difference photocurrent is then: 

h+i - h-i = 2 (^ajai^o + SiaLo) cos (fit). (S4) 

The cos (fit) modulation can be identified and de¬ 
modulated. By varying the phase of ano one can access 

{a\ ± hi) from the average values and (^a|±ai^ ) 
from the observed variation. From these quantities the 
covariance matrix can be reconstructed as in Ref. H. 



FIG. SI: Experimental scheme for measurement of co- 
variance matrix, a) Using phase delays {4>i, 4>2) and in¬ 
terfering on a beam-splitter generates the helds a± = 
± a 2 e*'^^)/v^. Different combinations of phases 
then allow access to the fields at as defined in the 
text, b) Interfering hi with a local oscillator of the form 
olo (e*^* -|- e“*^*) allows generation of the helds a±i = 
[hi ± olo (e“^* -I- e“*c^*)] /v^. The differences of the corre¬ 
sponding photocurrents gives access to the quadratures of hi. 

2. Analytic solution of the first and second order 
correlator evolution 

To calculate the fidelity of polaritonic CNOT gate as a 
function of decay rate F and pure dephasing rate F p , 
we derive the analytic solutions for amplitude and phase 
expectation values, as wells as higher order correlators. 
The solution of Eqs. (4)-(7) in the main text reads: 

(gi(t)) = (gi(0))e(^-i’-r-)*/(2?^), (S5) 

m)) = (^(92(0)) + ^(9i(0))f) , 

(56) 

(pi(t)) = (^(p^(o)) _ ^(p2(0))t) , 

(57) 

(P2(t)) = (?52(0))e(^-r-r-)‘/(2'^), (S8) 

where ((?„(0)) and (pra(O)) represent the initial amplitude 
and phase mean-field values. 

The evolution of second order correlators can be de¬ 
rived from Eq. (3) as: 

^ (1 “ 2 (^ 1 )) + — (1 + ‘2.{q\)) (S9) 

-*rp {{ql) - (S>1)) , 

= 2^J(g2) + i(p _ r - 2rp)(<jig2), (SIO) 

= 4*J(gig2) + f (1 - 2(92^)) (Sll) 

+ Y + ^^^2)) - *rp {{ql) - (P2)), 







+ f (1 - nfi)) (S12) 

+ ^(l + 2(p;))-ff, ((#?)-(©), 

+ i{p-r- 2r,)(p,re), (Si3) 

= f (l-2{®) + ^{l + 2(p=)) (S14) 

- *rp {{pD - {q^)) , 


= i{P-r-2Tp){qiP2), (S15) 

= -2iJ{qip2) + — —(2(giPi) - i ), 

(516) 

= 2iJ {{qipi) - {q 2 P 2 )) + i {P - T - 2rp) {q 2 Pi), 

(517) 

= 2jJ(giP2) + — —( 2 ( 92 ^ 2 ) - i) ■ 

(518) 

Eq. (IS^ - (IS18I) can be solved analytically for Ep = 0 
and r P, giving: 


dt 

dt 


j6<™> 

dt 


dt 


(flit)) 

{qi{t)q2{t)) 


mt)) 


ci + ((g?(0))-ci)e(^-r)*A (S19) 

C2 + ((171(0)92(0)) - C2) 

+ (c3 + ^(9?(0))) (S20) 

, , 4J((9l92(0)) - C2)^ 

C 4 + I (92 (0)) ~ C 4 H- - - 1 

2 j[c 3 h+ 2 J{qf{ 0 ))) „\ 

+ h? 

(S21) 


{m) 


m)) 


Cl + ((75i(0))-Cl)e(^-I’)‘/^ (S22) 

-C 2 + ((75i ( 0 ) 712 ( 0 ))+C 2 )e(^-n‘/^ 

- (c3 + ^(P^(0))) (S23) 

, „ 4J((piP2(0)) +C 2 )^ 

C4 + I (Pl(0)) - C4- - - 1 

2j(c3?i + 2J(p2(0))) 2 ^ (p-r)t/h 

+ ;i 2 * j e 


{qiP 2 {t)) = (9iP2(0))e(^ 

(9iPi(0) = 


(S25) 


+ (9iPi(0) - o - 


^ _ 2J(9iP 2(0))^ ^ (p_r)t/?i 




(S26) 


(<72P2(i)) = 


( 92 ^ 2 ( 0 ) 


_ 2J(9iP 2(0))^ ^ p(p-r)t/?i 

(S27) 


2 ;i 

2J 




{q2Pi{t)) = ( (92Pi(0)) + ((9ipi(0)) - ( 92 ^ 2 ( 0 ))) t 


4J^(9iP2(0))^^ ^ 


where the coefhcients are: 


Cl = 


C2 = 


C 3 = -2ci 

C 4 = Cl + 


p + r 
2 (r-p)’ 

2ci J 

r- p’ 

j 
h' 

4c2J 


r-p 

Separate equations hold for the special case E = P. 


(528) 

(529) 

(530) 

(531) 

(532) 


3. Multi gate realization: detailed description 

In the main text of the letter we described the scheme 
for a multi-CNOT gate, which can operate on different 
pairs of continuous variable modes, encoded by the pla¬ 
nar wavevector. The following can be done for instance 
in the spinfor dipolaritonic system, where direct exciton 
(DX), indirect exciton (IX) and cavity photon (C) form 
three distinct dipolaritonic modes UP, MP, and LP, with 
associated parametric scattering between them. Alter¬ 
natively, one can envisage a similar system, where two 
different direct exciton modes DXi and DX 2 are coupled 
to the same cavity mode C. 

We start the description from a generic dipolaritonic 
Hamiltonian [^, 

77dip=Po + Pi„t, (S33) 

where we separated the linear Hamiltonian of bare modes 
and associated couplings, Pqj and interaction Hamilto¬ 
nian Pint coming from exciton-exciton interactions. 

The linear part reads 

Po = X! ^ PDX,k,7^kjfek,7 (S34) 

k,7=± k.7=± 

+ X! -^IX.k,i4jCkj + ^C-Dx(dk,i4j + 4j«k,7) 

k,j = ± kj = ± 

- X! '4DX-Ix(akj4,4 + 4,i®kj), 
k,j = ± 


(S24) 
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where b]i,j (blj), Ckj (clj) correspond to 

annihilation (creation) operators of the cavity photon, 
direct exciton, and indirect exciton modes with k planar 
wave vector and j = ± circular polarization, respectively. 
The dispersions of the C, DX, and IX modes are Ec,k,j, 
EuxM.j, and ifix.kj', respectively, and we consider cir¬ 
cular modes to be degenerate in energy, i5.,k,-i- = £’-.k,-- 
f^C-DX corresponds to the exciton-photon coupling con¬ 
stant, and Jdx-ix denotes tunneling coupling between 
adjacent quantum wells. 

The interaction Hamiltonian can be written as 

"Hint = ^dx-dx + "Hix-ix + 'Hdx-ix, (S35) 

where we separate contributions emerging from direct ex¬ 
citon interaction, indirect exciton interaction, and direct- 
indirect cross-Kerr interaction. Starting the from conven¬ 
tional DX-DX interaction, the Hamiltonian reads 

■Hdx-dx = '^?°('^’^'’'^)^k-q,+^k'-rq.+^k,+^k',+ 

k.k',q 

+ «?°(k,k',q)6j,_q_6j,,+q_6k,-Sk',- (S36) 

k.k^.q 

+ '^?°(k.k',q)6j,„q^+6k,+&k-+q,_&k',- 

k.k^.q 

+ '^?°(k.k',q)6j,_q^„Sk,-&k.+q^+&k'.+ , 

k.k\q 

where corresponds to the triplet or co-circular 

Coulomb interaction between direct excitons of the same 
spin, and is the singlet or cross-circular interaction 
between direct excitons of opposite polarization. We note 
that for relevant momenta of the polaritonic system the 
exciton-exciton interaction constant is momentum inde¬ 
pendent and can be typically treated as a constant. 

In the following we are interested in particular wave 
vectors of the particles, namely the one where scattering 
happens between initial states of k 2 and ki, and final 
states k'l and k '2 chosen according to momentum con¬ 
servation. Using this labelling, Eq. (IS36I) yields 

'HdX-DX = Q^?°^k^, + ^k;. + ^k 2 ,-l-^ki.+ 

+ ,-^lc;,-^k2,-6ki.- (S37) 

+ “?°^k^,+^k2,-efck'^„&ki,- 

+ _Sk2,-6l,,+Sk,.+ . 

The interaction terms for indirect excitons can be written 
in a similar fashion, giving 

Hix-ix = Q!l^c^.+Ck'^,+Ck 2 ,-i-Cki,-i- 

+ “”4^.-4i,-Ck2.-4i,- (S38) 

-|- a2 _Ck2,-Ck'^__|_Cki,+ ) 


with triplet/singlet interaction between spinor indi¬ 

rect excitons. The cross-interaction between direct and 
indirect excitons reads 

HoX-lX = Q^?^^k^, + Ck'^, + ^k2. + Cki,-|- 

+ -^ki-^k2.-4i,- (S39) 

E OI 2 ^k^, + ^k2, + Ck/^__Cki,- 
T 0^2 ^k2, —^k2, —^k^,-t-^ki, + ; 

where of 2 corresponds to direct-indirect exciton 
Coulomb interaction for same and opposite spin projec¬ 
tions ±1. 

Next, we perform the transformation of the linear 
Hamiltonian (IS34p to the diagonal basis of dipolariton 
operators Aj {j = 1,2,3): 


b-k,± = Vii^k,±Ai^k,± + V2l,fc,±^2,fe,± + V3i^k,±A3^k,±, 

(540) 

bk,± = Vi2^k,±Ai^k,± + V22,fc,±242,fc,± -I- V32,k,±A3^k,±, 

(541) 

Cfc,± = Vi3^k,±Ai^k,± + V23.fe.±^2,fe,± + V33_fc_±243_fc_±, 

(542) 

where Vij are the matrix elements of eigenvectors, cor¬ 
responding to dipolariton Hopfield coefficients. Without 
limiting the generality we consider Vij to be real for any 
i and j. 

The transformed linear part reads 


Hq — J, j_24j,k,±. (S43) 

i,k.± 


The eigenenergies possess rotation symmetry, i?j(k) = 
Ej{k), and thus the Hopfield coefficients are spin inde¬ 
pendent, Uij,k,± = Vij, where we omit the momentum 
index for brevity. 

We proceed with the transformation of the interaction 
Hamiltonian by successive transformation of its parts. 
For instance, the direct exciton interaction part can be 
rewritten as 
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T^DX-DX = .+ + V22AI^,,^ + . + )(Fi2il,^,+ + + ^^32^^,,+ ) (S44) 

(^12^1,k2,+ + 4^22^2,k2,+ + 4^2^3,k2, + )(4l2^1,ki,+ + V22^2,ki,+ + 4^32^3,ki.+ ) + (4l2^| _ 

+ 422^2,k^,- “4 4^32^3 i(;^,-)(4l2^1^k'j,- "4 ^2^2,k^.- “4 432^3 k'^ _-) (4l22il,k2,- + 4^22^2,k2,- + 4^2213,k 2 ,-) 

(4l22ll,ki,- + 442^2,ki,- + V32^3,ki,-) + (4l2^i,k^,+ ”4 ^ 22 ^ 2 ,k^,+ ”4 ^ 2^3 k^+) 

(4l2^1,k2,+ + 4^22^2,k2,+ + 4^2^3,k2, + )(44l2^J k'j.- “4 ^2^2,kj,- “4 4^2^3 k'^^,-) (44L2^1.ki,- + 4^222^2,ki,- + 4^2^3,ki,-) 
+ Ck°'^(Vi2^3 k^ - + 4^22^2,k^,- ”4 ^32^3 k^ -)(4l2^1,k2,- + 4^22^2,k2,- + 4^32^3,k2,-) 

(Vi 2 ^J k'j^,+ + ^ 22 ^ 2 ,ki,+ “4 ^ 2^3 k'j__|_)(4^12^1,ki,+ + 4^22^2,ki,+ + 4^32^3,ki,+)■ 

I- 


To proceed, we recall the relevant quantum and classical 
modes required for the operation of a serial multimode 
CNOT gate shown for example in Fig. 3, main text. As¬ 
sociating the UP, MP, and LP modes with Ai, A 2 , and A 3 
modes, the generic interaction Hamiltonian (IS44I) can be 
accommodated to our needs using the set of relabellings: 

Ai,ki,-i- '0up,-i-i 2ii_ki,- 0, (S45) 

Al^k2,-|- Oj 2ll_k2,- 0) 

A2.ki,+ Si,-1-1 2l2,ki,- ^1,-! 

A2,k2.+ ^ ^2,+ , A2,k2,- ^ S 2 ,-, 

A3,ki,+ 0, As^ki,- '0LP,-i 

A3,k2.+ '-t V’LP.+ i A3_k2.- '-t V'LP.-i 


and again the final state (primed) wavevectors are chosen 
according to momentum conservation rules. These imply 
that ki -f k 2 = k 2 -f k'^ for parametric coupling processes 
and k 2 —ki = k^ —k'^ for linear coupling processes, which 
can always be satisfied by choosing wave vectors for the 
classical modes appropriately. 


Using Eq. (IS45I1 we can rewrite Eq. (IS44() as 


iy22a\ j^ + 4^32'0LP,-|-)(4l2'0UP,-|- + 4422S| _|_) (4^2 V'UP.-I- + 422201,+) (V 22 02,+ + 4(32V'LP, + ) 
+ 4232'0lP,-)(4422O3 _ + 4232'0L*p, _ ) (1222 Ol,_ -|- 4232'0 lP,-) (4^2202,- + l232l/'LP,-) 

+ 4232V'lP,+ )(44i2V'UP,-|- + 422201,+ )(422A| _ -|- 1232i/’lP,-)(422A 2 ,- + V^2i^'LV,-) 
LP,-)(42220i,_ -f 1232 V'lP,-)(4A2'*/’uP,+ + 4220^ _|_)(422a2,+ + 4232V'LP,+ )- 

I- 


(S46) 


^DX-DX — Oi 

+ a?°(U22a^,_ 

-ba°°(U22a^_+ 

+ C1?°(422202,_ + 1232 V' 

Taking only energy conserving terms, we arrive to the 
Hamiltonian derived from the DX-DX interaction: 

’^DX-DX = a?°(442|lA2432V'UP,-i-V'LP, + a 2 ,_,_a| -I- h.c.) 

(S47) 

+ a?°( 4 ^ 22432 ' 0 LP,-' 0 LP,-O^__ai,_ -b h.c.) 

+ a°°( 4 " 224 ^i 24 A 2 i/'up.-iV'LP,-aJ_+a| _ -b h.c.) 

+ Q;^^(4222432^LP,-'0LP,-Oi,+a2,+ -b h.c.) 

+ a°°( 4 " 22442 ^LP,+V'LP,-al _ai,+ -b h.c.). 

A similar procedure can be applied to IX-IX and DX- 
IX interactions, leading to the transformed interaction 


Hamiltonian of the form: 


■^int = (Ai/'up,-e'0LP,+4,+4,+ + ^-c-) (S48) 

-b (Hi/>LP,-V'LP,-a 2 ,_ai,_ -b h.c.) 

+ (C'' 0 UP,-i-V'LP,-a 2 ,+a| _ -b h.c.) 

+ (£>i/>lp,-V'lp,-4,+«i,-i- + ^-c-) 

-b (Ei/)£p,+' 0 LP,-ol__ai,+ -b h.c.), 
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where we defined the constants 


hand side of Eq. (IS56I) : 


A = a^^V22^12V32 + a^^l23Vl3V33 + a^^^V22V23 V 12 E 33 , 

B = + al^lVsl + (S49) 

C = Q;^^E 22 ^12^32 + 02 ^^ 23 ^ 131/33 + 0 ^^ 22 ^ 23 ^ 121 / 32 , 

7 -) _ ^DDt/2 J/-2 I ^IIt/'2 t/2 I ^DIt/'2 t/2 

11 — 02 E 22 V32 + 02 V23 E 33 + 02 *123 *421 

E = + al^lVsl + a^^V22y2-iV33- 


To achieve the desired interaction between quantum 
modes we perform a polarization transformation from 
the circular to linear basis, accounting for the TE-TM 
splitting, which is present in exciton-polaritonic systems. 
The transformation reads: 


1 


a 

a 


V2 

1 

’■"71 




(fin.TMe — ifin.TEe , 


(S50) 


where (f>n (n = 1,2) represent angles encoding momen¬ 
tum states. We are mainly interested in the TM inter¬ 
action terms (target quantum modes oi^tm and 02,tm), 
while TE modes as well as cross-terms can be disregarded 
in the case of large TE-TM splitting (see the discussion 
in the next section). Also, we write the classically driven 
modes ijjj explicitly as complex numbers with absolute 
value of and phase (j)j^ namely: 

V'up.+= ^up,+e*‘^^^-+, 4’UP.+ = «'LP,+e*‘^^^'+, (S51) 

V'LP.- = 4'LP,-e*'^*i^’-, V'LP.- = «'LP,-e*‘^^^’- • 

Performing the transformation of Hamiltonian (IS48I) , we 
get: 


"^int = 

C 


A 




-'I'up,+'I'Lp,-e 


*(0UP, + -|-0LP,- —202-|-2<^i) 


®2,TM®1,TM 


-b 


B 


d^LP -'f'LP _e*(''^*'*"'-”7p,-+202-2b'i) 


B , 


d^LP -'I'lP _e®(‘^^’"'““’^LP,--202-l-2<^.i) 


^2,TM^l5TM 




b(0UP, + -|-b'LP, + -2</'2-20l) _ J 

(S53) 

b('/'UP, + -l-</’LP,-— 2<^2-|-20i) ■ 

— t, 

(S54) 

b(0LP,--?!>LP,-+202-2<^i) _ J 

(S55) 

b(0LP,--b>LP,--202-1-201) _ 

(S56) 


The first system of equations (with plus sign) can be 
satisfied for 


<p 2 = (l>i + n 7 r/ 2 , (S57) 

4up.+ = —0LP,- + (2n -b l) 7 r/ 2 , 

4>hp,+ = 4>hp,- + 4411 

0LP.- = 4>PP- - (2n -b l)7r/2, 

and we have the freedom in choosing 4 lp,- phase. For 
this choice of phases, the interaction constant shall be 
tuned to (A4 'up,+'I'lp,-i- + C'4^up,+4'lp,-)/ 2 = (H -b 
-D)<p,_^lp.-/2 = J. 

The second system of equations (with minus sign) can 
be satidfied for 

<(>2 = <('1 + ?T- 7 r/ 2 -b 7 r/ 4 , (S 58 ) 

(()up,+ = -4>PP- + (/I + l)7r, 

4lP.+ = 0LP,- + 4(/)i, 

<(’lp,— “ 4 ^lp,— ~ j'ii 

with interaction constant tuned to (A4'up,+4'lp,-(- + 
C^up,+4'lp,-)/2 = iB- i^)<p__^LP.-/2 = J. Ac¬ 
counting for the possibility to modify the system on the 
fly with adjustable pumps, we can in principle organize 
a sequence of gates between eight different momentum 
modes $1 = {4i,</>i + 7 *‘/ 4 ,..., </>! -b 77 r/ 4 } characterized 
by wave vectors ki. Finally, coupling to another sub¬ 
space of continuous wave modes $3 defined by ip'i can be 
done with lower fidelity and exploiting error correction 
afterwards. 

Given the versatility of the system, we can efficiently 
control the couplings by classical drive amplitudes, and 
thus arrange the H = 2 Jgi_TMP 2 .TM Hamiltonian for 
original interaction constants a"™ {nm = DD,H, DI), 
including both triplet and singlet interactions. Addition¬ 
ally, the phase conditions set the energy shift 

to E4 'lp,-i-'I'lp,-, implying that it shall be minimized 
by weak drive conditions for V'lp,-!- and V'lp.- classical 
modes, and stronger pumping of 4 up,+ and '^pp - modes. 


where the first two terms correspond to useful parametric 
and linear coupling terms, and the third term represents 
an additional energy shift of one of the modes. Finally, 
to reduce the system to the required qiP 2 type of interac¬ 
tion, the phases of the classical drives can be adjusted to 
make each term in (IS52I) purely imaginary. This can be 
satisfied with two sets of conditions defined by the right 


4. Multi gate realization: TE mode influence 

Now, let us return to the question of the rotating 
wave approximation (RWA) validity for the TE-TM 
cross-interaction terms. So far we kept only TM modes, 
assuming that the TE-TM interaction is strongly sup¬ 
pressed due to the energy shift of TE modes, i.e. TE-TM 
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FIG. S2: Influence of a spurious TE-TM couplings on the 
fidelity of multimode CNOT gate. The minimal fidelity is 
shown as function of additional coupling rj for various mode 
occupations {{x)'^ = 1 to {x)^ = 10), and three values of 
dimensionless TE-TM splitting. 

splitting. However, if the splitting Ate-tm is small com¬ 
paring to other relevant energy scales (interaction con¬ 
stants and decay rate), this assumption becomes invalid. 
To test the limits in which auxiliary TE modes can be 
neglected, we recall that together with useful terms ap¬ 
pearing in (IS52I) . various spurious terms appear, namely: 

e22d2,TE®2.TE, (ei2apTE®2,TE+h-C.), 
(b21®l.TM®2,TE~l”^'*")’ (^120^ rpjyjClj rpj^-|-h.c), 

(C 2 ia 2 TM^i.TE+h.c), etc. Here constants tij, rjij 
and Cij denote generic coupling represented by functions 
of bare interactions, Hopfield coefficients, and classical 
drive amplitudes. The first three terms act fully in the 
extra subspace of TE modes and are irrelevant for our 
considerations. However, the terms of fourth, fifth and 
so on type produce the parasitic rotation for the system 
modes ai/ 2 ,TMj or alternatively effective leakage to 
the additional TE mode subspace. The related fidelity 
degradation of the gate then depends on the values of 
couplings, mode detuning Ate-tm, and decay of the 
mode. 

To quantify the fidelity change due to additional cou¬ 
plings, we refrain from considering a particular system 
with defined coupling, but characterize the generic influ¬ 
ence of e.g. rj type of the coupling. For this, we consider 
Hamiltonian: 

T-L = J(— mi,TMa2,TM + rpjy; — rpj^a2,TM 

(S59) 

fl- ^2 rpj^Ui^TM) “t” ^(^2.TM^1.TE h.C.') Au| -pE; 


where A is a generic TE mode shift defined by TE- 
TM splitting and nonlinear contribution. Here, for the 
sake of simplicity let us rename modes as oi.tm = di, 
02,TM = 02, and ui^te = 03. Next, Eq. (IS59I) can be 
rewritten using position and momentum operators qj and 
Pj associated to each mode, which gives: 


n =2Jqip2 - r]{q2P3 + P2q3) + ^{ql +pI)- (S60) 


Deriving the equations of motion for average amplitudes 
{qj), (pj) {j = 1,2,3) and associated correlators, we can 
calculate the fidelity for CNOT gate acting in the {oi, 02 } 
mode subspace as a function of dimensionless parameters 
p/J and A/J. 


The results are shown in Fig. IS2[ where we consid¬ 
ered the cavity decay rate to be small, T/ J <C 1. We see 
that even for spurious interaction constants p being com¬ 
parable to coupling J, the degradation of fidelity can be 
suppressed by the shift of TE mode. In particular, taking 
J = 0.1 meV for a dipolaritonic system, and assuming 
realistic 0.5 meV TE-TM splitting one can achieve the 
0.99 fidelity commensurable with previously anticipated 
degradation due to decay of the cavity mode. 
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